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We use the radial gauge to calculate the recently proposed 
ansatz for the physical electron propagator in such effective 
models of strongly correlated electron systems as the QED3 
theory of the pseudogap phase of the cuprates. The results 
of our analysis help to settle the recent dispute about the 
sign and the magnitude of the anomalous dimension which 
characterizes the gauge invariant amplitude in question and 
set the stage for computing other, more physically relevant, 
ones. 



I. INTRODUCTION 

As a generic property, the one-dimensional Fermi sys- 
tems with short-range (screened) repulsive interactions 
routinely demonstrate algebraic decay of all the cor- 
relation functions governed by non-universal (coupling- 
dependent) anomalous exponents. 

A possibility of the emergence of a similar behav- 
ior, commonly referred to as the "Luttinger liquid", in 
higher dimensional strongly correlated electron systems 
has been extensively discussed in the recent literature. 

Thus far, however, no solid consensus has been reached 
even on the necessary criteria that have to be fulfilled for 
the Luttinger behavior to set in, much less on whether 
or not it occurs in any specific example of a strongly 
correlated electron system. It was largely for this rea- 
son that the attention has recently been drawn to the 
class of effective models described by (possibly, spatially 
anisotropic and/or Lorentz non-invariant) deformations 
of the standard action of Quantum Electrodynamics. 

Motivated by the puzzling properties of the quasi- 
two-dimensional high-temperature copper-oxide super- 
conductors, most of the interest has been focused on 
the three-dimensional case described by either the or- 
dinary (parity-even) QED3 or the abelian 3D Chern- 
Simons theory 1,2 , where the finite density problem of 
non-relativistic (massive) fermions has become the main 
subject of the scrutiny. However, the latter was found 
to fall into a rather different class of non-Fermi liquids 
which bear little resemblance to the ID Luttinger liquid 3 . 



Recently, the idea of the conjectured Luttinger-like 
behavior has been rekindled in the recent theories of 
the pseudogap phase of the underdoped superconducting 
cuprates 4-8 . Albeit describing rather different physics, 
all these approaches resort to the same effective descrip- 
tion in terms of the pseudo-relativistic QED3 theory of 
the gapless nodal fermion excitations which retain their 
d-wave symmetrical gap well above the critical temper- 
ature T c regarded merely as the onset of global phase 
coherence. 

Above T c the fermionic excitations experience strong 
scattering by both thermal and quantum fluctuations of 
an incipient ordering, such as a flux of the gauge field 
measuring a local spin chirality 4,5 or vortex-antivortex 
pairs of the pairing order parameter 6-8 . The latter sce- 
nario has recently received a new experimental support 
from the observation that the vortex matter is present 
at temperatures well in excess of T c , as revealed by the 
measurements of the Nernst effect 9 . 

As another important ingredient, the QED3 theory of 
the pseudogap phase was aimed at explaining the ubiq- 
uitous destruction of the coherent quasiparticles above 
T c which was observed in angular-resolved photoemission 
and tunneling experiments. 

To this end, the authors of Ref. 5 conjectured that the 
electron propagator in question may, in fact, exhibit the 
Luttinger behavior 



G(x) oc i/|2:| 3+ ' ? , x = 7 M : 



(1) 



characterized by a positive anomalous exponent r\ > 0, 
and they also attempted to fit the ARPES data, while 
claiming a good agreement with experiment (unless ex- 
plicitly stated otherwise, throughout this paper we use 
the notation n = "f^n 
the Dirac matrices 7^). 

The conclusions drawn in Ref. 5 were based on the use 
of the following heuristic form of the gauge-invariant elec- 
tron propagator 



M for vectors contracted with 



G(x — y) = (Q\ip(x) exp 



i / dz^A^z) 



mm, (2) 
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where the line integral was taken along the contour T 
chosen as the straight-line segment connecting the end 
points. 

Later on, the calculations of Ref. 5 were carried out 
by a number of other authors, and the results for the 
anomalous exponent appeared to vary not only between 
the different authors {rj = 32/3n 2 N 5 vs -32/3tt 2 N 10 ) 
but also from one to another work of the same authors 
{r} = -16/3ir 2 N 6 vs 16/3Tr 2 N n and r\ = 32/3tt 2 N 7 vs 
-64/37r 2 iV 12 ). 

While some of the calculations were performed in the 
conventional covariant gauges 10 ' 12 , other authors made 
use of the potentially problematic axial gauge ((x — 
y) fJ _A fJ _(z) = where x,y are arbitrarily chosen points 
which are taken to coincide with the end points of the 
contour r) 5 ' 6 ' 10 which spurred a debate over the issue 
of a true (vs limited, see 13 ) gauge invariance of Eq. (2), 
as opposed to its surrogate functions proposed in 11 (see 
Summary for a more extended discussion). 

While seemingly being an issue of secondary impor- 
tance, a proper construction of the physical electron 
propagator is, in fact, imperative, as far as ascertain- 
ing the status of the conjectured Luttinger-like behavior 
in the QED-like theories is concerned. 

In light of the present controversy, in this paper we un- 
dertake yet another attempt to settle the dispute about 
the physically motivated form of the electron propagator 
and the actual value of 77 (if any) by resorting to the so- 
called radial (Fock-Schwinger) gauge ((z — x)^A^(z) = 
and x is an arbitrary fixed point). The radial gauge is 
known to be free of the potential problems that might ex- 
ist in the axial gauge, which, according to some authors, 
may even require one to introduce ghost fields 14 . In ad- 
dition, we also set out to explore the dependence of the 
previously conjectured form of the electron propagator 
(2) on the choice of the contour T. 



II. GAUGE INVARIANT FERMION 
PROPAGATOR IN THE FOCK-SCHWINGER 
GAUGE 

We start with the 3D relativistic theory of massless 
Dirac spinors coupled to a massless U(l) gauge field, 
whose Euclidean action is 



N 



(3) 



N = 2. Nevertheless, in what follows we choose to treat 
N as a parameter that can assume arbitrary values, de- 
pending on the problem in question. 

The dynamics of the U(l) gauge field is generated 
by the effective action obtained after tracing out the 
fermionic degrees of freedom 



S ctf [A] = IjcTxJ d n yA,{x)D-l(x-y)A v 



(y) 



(4) 



where the dots stand for the higher order terms (non- 
Gaussian) which we hereafter neglectm as is done in all 
of the previous works on the subject. Intending to sub- 
sequently use the method of dimensional rcgularization 
for evaluating Fcynman diagrams we formulated Eq.4 in 
n = 3 — e dimensions. Also, in Eq.4 we neglected the 
bare Maxwell term ~ {d^A v — d u A^) 2 which turns out 
to be irrelevant in the low-energy, long-distance limit. 

The previously proposed candidate for the physical 
electron propagator (2) studied in 5-7, 10-12 can be cast 
in the following form 



G inv (x-y) = {G[x, y; A] exp 



dz» A^z) 



), (5) 



where G[x, y; A] is a fermion propagator for a given fixed 
configuration of the gauge field A^, and the brackets 
stand for the (normalized) functional average over the 
gauge field which is described by the action Eq.(4). 

In the Euclidean momentum space, the kernel D~* of 
the quadratic operator has the form 



(6) 



n-l/ \ Ny/q 2 I q^q v 
(q) = — ^— 



Introducing a source field 



1 

J»{z) = {x-yy j da6 n (z-y-(x-y)a), 



(7) 



we can write the straight-line integral which appears in 
Eq. (2) as 



X 

J dz^A^z) = J d n 



zJ»(z)A»{z). 



(8) 



where -0 = V^7o and the iV-fiavored Dirac fermions are 
described by four-component bi-spinors which belong to a 
reducible representation of the 7-matrices satisfying the 
Clifford algebra {7^,7^} = 25^ = 0,1,2). The 

latter can be chosen in the form of the direct product 
7^ = 03 (g) (03, 02, <ti) of the standard triplet of the Pauli 
matrices a^. 

In all of the above mentioned condensed matter- 
inspired QED-like models the number of fermion flavors 



In the Fock-Schwinger gauge 

(x - x )^A f _ l (x) = 



(9) 



the line integral in Eq. (5) vanishes if one chooses the 
reference point x at the "center of mass" x = X = 
(x + y)/2 (for the proof, see, Appendix A which also 
contains a derivation of the photon propagator in the FS 
gauge D™( Zl +X, z 2 + X)). 
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We compute the first order 1/N correction to the 
fcrmion propagator (3) by expanding the inverse Dirac 
operator l/(d-iA) = l/d+ (l/<9)(ii)(l/d) + . . . in pow- 
ers of A^(z), thus obtaining 



G inv (x, y) = S(x -y) + SG(x, y), 



(10) 



where S(x — y) and 5G(x, y) are the bare fermion propa- 
gator and the correction to it, respectively. The latter is 
given by the expression 

SG{x, y) = - J d n Zl J d n z 2 S(x -z x - X)Y 

xS( Zl - z 2 ) 1 v S{z 2 + X- v)D™( Zl + X,z 2 + X), (11) 

where it is natural to decompose the propagator for 
the gauge fields in the Fock-Schwinger gauge into four 
contributions, 



D™(x,y)=J2D$(x,y), 



(12) 



i=i 



as is done explicitly in Eqs. (A20) to (A23). By choosing 
the reference point in the Fock-Schwinger gauge to be the 
center of mass X, one verifies by inspection of Eqs. (A20) 
to (A23) that the contribution SG can only depend on the 
relative coordinate x — x — y: 

SG(x) = - J d n z x J d n z 2 S(x/2 - Zl yf 

xS( Zl - z 2 )YS(z 2 + x/2)D™(z x + X,z 2 + X). (13) 

Taking the Fourier transform of Eq. (13) FT[5G(x)] = 
i8G(p), FT[S(x)} = iS(p), S(p) = 1/p, we obtain 

xS(p - qi /2 + q 2 /2)YS(p - qt/2 - q 2 /2)Df v {q u q 2 ). 

(14) 

In the above equation, the index i = 1,...,4 labels 
the Fourier transforms D^l(qi, q 2 ) of the components 
D$(Z! +X,z 2 + X) given by Eqs.(A20)-(A23) 

£>#(«i,<fe) = (2n) n 6 n (q 1 +q 2 )6 tlv D(q 1 ), (15) 



O i l rJ(<h--'h)= I daq^-^-D^(q u q 2 ), (16) 



D$(qi,<h)= J ' d(3q 2v -^-D^(q 2 , qi ), (17) 





1 1 

D ^KlU I2) = J da J d(3q lfl q 2v 



_d d_ 

dqix dq 2X 



D^\q u q 2 ), 



(18) 



where 



D(q) = 



N 

D {a) {qi,q 2 ) = (2n) n S n ( qi +aq 2 )D( qi /a), 
D^\q u q 2 ) = (2w) n 6 n (l3q 1 +aq 2 )D(q 1 /a). 



(19) 

(20) 
(21) 



The first contribution to the anomalous dimension stems 
from SG^^p) given by Eqs. (14) and (15). Since the 
Fourier transform of the propagator Eq.(l) is propor- 
tional to p/p 2 ^ 11 it proves convenient to multiply SG^ (p) 
by p in order to deduce the exponent rj which we are after. 
Now taking the trace we arrive at the following result 



i T r 

4 



SG (1 \p) 



d n q 



(27T 



-Tr[pS(p)rS( P -q)r 



xS(p)]D(q) = - 
3tt 2 N e ' 



8(n- 2) 
N 



f d n q 



p-q 



/ tp 2 (p — q) 2 



(22) 



where we used dimensional regularization near n = 3 — e 
in the divergent integral 



p' - p-q 



d n q 

Wr y/?p 2 (p-q) 2 



J_\p\_ 

6ir 2 e 



(23) 



Notably, Eq. (22) coincides with the result for the anoma- 
lous dimension of the ordinary (gauge-variant) fermion 
propagator performed in the covariant Feynman gauge. 

Next, we compute the term 5G^ 2 '(p) given by Eqs. (14) 
and (16)) 

«' 2 >( P ) = /«o./^L05( P + „/2 + E /2) 



x^S(p - Ql /2 + q 2 /2) 1 v S{p - Ql /2 - q 2 /2) 

xq lfl ^-(2n) n 6 n ( qi + aq 2 )D( qi /a). (24) 
oqiv 

Integrating by parts we cast SG^ (p) in the form 
l 

SG^(p) = -JdaJ 00(^1^) 



x^)- 8 



a dq lv 

qi 



S(p 

12, 



(25) 



where we made use of the Ward-Takahashi identity 
(WTI) for the bare propagators 
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S(k + q)qS(k) = S(k) - S(k + q). (26) 
Multiplying Eq. (25) by p and taking the trace we obtain 



1, 



f>5G {2) (p) 



da 



d n qi d n q 2 



(2?r) n (2tt)" 
u 

x(2TT) n d n (q 1 +q 2 a)D(q 1 /a) 

(n-2) (2p 2 -f>- 9l ) 



x< - 



+ 



{ 2 (p-gi/2 + g2 /2) 2 (p-<Zi/2- g2 /2)2 

(n-2) (P-91+P-52) 



2 (p-g 1 /2- <Z2 /2) 2 (p + gi /2 + g2 /2) 2 



• (27) 



Being an odd function of the momenta qi^ the last term 
in Eq. (27) vanishes. Thus, we obtain 



-Tr 



p5G {2 \p) 



4(n-2)2" f f d n q 



N 



(27r) n 



(p 2 +p- qa) 



q 2 (p + q(l + a)) 2 (p - q(l - a)) 2 



(28) 



The momentum integral can be evaluated by virtue of 
the Fcynman parameterization 



h 



d n q 
(2tt)" 



(p 2 + P • ga) 



V(p + g (l + a)) 2 (p- (? (l-a)) 2 



(l-a 2 ) 3 4 



1 1 

:/*/ 



o o 



(1 — (2xy — x)a 



-(l-x)a 2 ) J 



dq 



„n-l 



2tt 2 [g 2 + cf' 2 



o 2 3 



(l-a 2 ) 3 4 



- dx dy 



-(l-x)a 2 



o o 

„(n-5)/2 



(1 — (2xy — a;)a 



6tt 2 ' 

where the argument c takes the form 



(29) 



2 

p X 



(I -a 2 ) 2 
+ (l-a;)a 2 ] . 



[l-x(l-2y) 2 + 2(1 -20(1-2^ 



(30) 



The leading divergence of the integral over a can be ex- 
tracted by making an approximation similar to that of 
Ref. 16 



„(n-5)/2 



(1 - a 2 ) 2 



(n-5)/2 



a; [1 - x(l - 2y) 2 + 2(1 - x)(l - 2y)a + (1 - x)a 2 } ' 

(31) 



In this way, we obtain 



h 



87r 2 (l-a 2 ) 



|P| 



(1-a 2 ) 



1(a), 



(32) 



where 



l l 
1(a) = j dx j dy 
o o 

(1 - x)-!/ 2 (l - (2xy - x)a - (1 - x)a 2 ) 
X [1 - z(l - 2y) 2 + 2(1 - x)(l - 2y)a + (1 - x)a 2 



(33) 



The expression for <5G^ 2 ' now reads 

l 



P 5G {2) (p) 



1(1) da 



4|p|" 



r 2 7V 



da 



/(a) 



1 



(l-a 2 ) 1 ^ 
4 



2(1 -a) 



l-e 



TVtt 2 



(34) 



where we used the integral 
HI) 



/ dy ~ttt 77T = 2. 







2[l + x(y-l)] 



(35) 



It can be readily shown that Eq. (34) is also identical to 
the result for Tr\p5G^ (p)}, given by Eqs. (14) and (17). 

Lastly, the expression for SG^ 4 \p) given by Eqs. (14) 
and (18) reads as 



f-^L [dr 



a 



n-2 



j- d n qi d n q 2 
J (2tt)" (2tt)" 



x^-l + f^-f-f) 

(36) 

where t = (3 /a. First, we make use of the WTI given 
by Eq. (26) to simplify the product of the free fermion 
propagators 

s( P + q i + f)^-| + f)^-|-f) 

= ^+| + f)-^-| + f) 



As a result, Eq. (36) assumes the following form 



(37) 



4 



1/0 





x(2^)"«5"( gi r + g2 )^( gi ) 



_d d_ 



(38) 



where we have dropped the last term in Eq. (37) which 
vanishes upon angular integration. Next, we pull the 
derivatives to the front of the integral 

W4), 1 = H1 / <fa Yd f dnqi 

W A8 PX 8 PX J a"" 2 7 T J (2tt)» (2tt)» 



4 dp A dp A 



x (27r) n J n (giT + <fc)D(«i) [S(p + gi/2 + q 2 /2) 
+ S(p- qi /2 + q 2 /2)} 



(39) 



and carry out the momentum integration over q 2 fol- 
lowed by rcscaling of the remaining momentum variable 
qi which yields 



o o 



d n q 
(2n) n 



D(q) 



x [S( P + q(l - r)/2) + S(p - q(l + t)/2)} 

1 1/q 



2 n d d [ da f 

nw,wJ^I dT [{l - T) 



l-n 



+ (1 + r) 1 -"] 
~ ^ 2 iV e ' 



(p + g) 



(2tt)" v^(p + g)2 



(40) 



where we have invoked Eq. (23) to compute the integral 
(p + tf) P \p\~ e 



(27r)» v^(p + q)2 6^ 2 e 
and used the n — > 3 asymptotic 

1 l/a 



(41) 



(42) 



Notably, a singular (1/e) term which is present in the 
gauge propagator D$(qi, q 2 ) and the origin of which is 
discussed in the Appendix gets canceled in Eq. (40). 

Combining the four contributions (22), (34), and (40) 
we finally obtain 



Tr \pSG(p)} = ^Tr [p (sG^(p) + 28G {2) { P ) 



+ 



SG {4 \p) 



32 \p\~ 
'Zir 2 N e 



(43) 



which implies that in the momentum space Eq. (2) ac- 
quires the form 



G inv (p) = S(p) + SG(p) =^[1~V - In H)] , (44) 

thus allowing one to read off the anomalous dimension 

32 



n : 



3tt 2 7V' 



(45) 



This result corroborates the earlier calculations per- 
formed in the covariant and axial gauges in the frame- 
work of both the path-integral approach of Ref. 10 and 
the direct perturbative expansion of Ref. 13 . 



III. GAUGE INVARIANT PROPAGATOR WITH 
SEMI-INFINITE STRINGS 

In this section we set out to investigate the depen- 
dence of the amplitude (2) on the choice of the contour 
T. Specifically, we consider the contour consisting of two 
(anti)parallel semi-infinite strings attached to the end 
points, in which case the source current J M (z) in the line 
integral (8) is given by the expression 

oo 

J"(z) = n»Jda[8 n (z-x-na)±8 n (z-y±na)], (46) 
o 

where is a unit vector in a direction of strings. The 
upper (lower) signs in (46) correspond to the cases of 
parallel and antiparallcl strings, respectively 

In the former case, despite the fact that the strings do 
not form a closed contour, the corresponding amplitude 
remains gauge invariant, as long as all the infinitely re- 
mote points can be compactified into a single one. This 
customary assumption always holds in the perturbative 
sector of the gauge theory where all the fields vanish at 
infinity. 

The correction to the free fermion propagator has a 
form 

6G(x -y) = ~ J d n z 1 d n z 2 J^(z 1 )D flv {z l - z 2 )J v {z 2 ) 
xS(x-y) + j d n z 1 d n z 2 S(x-z 1 )^S(z 1 -y) 
xD^Z! - z 2 )J u (z 2 ) - J d n Zl d"z 2 S(x - z^-y^ 



x S(z! - z 2 )-y fJ ,S(z 2 - y)D liV (z 1 - z 2 ). 



(47) 



Making use of the amplitude G(x — y) being explicitly 
gauge invariant, we choose to compute it in the Feynman 
gauge where the gauge propagator takes a particularly 
simple form 



(x*y 



A = 



4 fi 3 - n r{^] 



N 



n±l 
7T 2 



(48) 
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The last term in Eq. (47) corresponds to the standard 
fcrmion self-energy in the Feynman gauge 

Ntt— (2n - 3) 3-n 

The first term in Eq. (47) containing two source currents 
is given by the integral 

5G^(x ~y) = -\ Jd n Zl d n z 2 J^D^Z! - z 2 )J u (z 2 ) 



= -\ Jdajd(3 







[{x-y + n{a- (3)f 



[{-x + y + n{a - ft)) 2 } — 



(50) 



The integral over (3 is convergent for 1 < n < 2. However, 
the remaining integration over a diverges with the upper 
limit L which we impose as a cut-off. 

Rescaling the integration variables a —* a\x — y\, (3 —* 
f3\x — y\ and introducing the angle 9 according to the 
relation cos 6* = (x — y) ■ n/\x — y\ we rewrite Eq. (50) as 



5G {1) (x-y) 



A\x - y 



3 — n 



L/\x-y\ oc 



da / dj3 



la-/?!"- 1 



[1 + 2cos6»(a -13) + {a- f3) 2 }^ 
1 

[1 - 2 cos 6{a - (3) + (a - /3) 2 ] V 



• (51) 



In the integral over f3, we consider separately the intervals 
from to a and from a to oo. First wc compute the 
integral from a to oo which, upon shifting the integration 
variable (3 — > (3 + a, acquires the form 



h 



dT 



1 



f 



(r 2 Y 



[1 - 2tcos(9 + t 2 ] i V 



(52) 



After exponentiating the denominators and carrying out 
the integral over r one obtains 



h 



r(^) 



dss 2 



*-2 



1 - e 



-ssin 



v^r(f -l) . 2 u _ a 

/ 2 (sm 2 ftV 2. 

r(s=i) v ; 



(53) 



In turn, the integral over f3 from to a takes the form 

a 

d(3 



[1 + 2cosft(a -(3) + {a- f3) 2 }^ 



d/3 



[l + 2a/3cos(9 + a 2 /? 2 ] 



(54) 



where we made a change of variables (3 — > a — /3, followed 
by rescaling /? — » a/3. 

Next, we represent the last expresion as the difference 
between the integrals taken from to oo and that from 
1 to oo 



dp 



[f3 2 + 2/3 cos ft + 1] — 



OO 

w 



oo 

-(a 2 + 2acosft' + l) 1 ^ J 



d(3 



[f3 2 + 2(3 cos 7 + 1] — 



r, (55) 



where 



cos 7 : 



a + cos ( 



\/a 2 + 2a cos 6* + 1 



(56) 



The integrals in Eq. (55) are evaluated with the use of 
the formula 



oo 

/ 



dx 



(x 2 ± 2xcos7+ 1)p 2 



= -(sin 2 7)2 ^ i 2^ 



r(p) 



(57) 



+ COS7F \ l,p; -;cos 7 J . 

Hence the sum of the two integrals reads 

7 2 (cos0) + / 2 (-cos0) = 
= (a + cosft)(a 2 + 2acos6> + l) 1 ^ 

(n — 1 3 \ 
1,— ^— ;-;cos 2 7 J +(cos6»^-cos6»), (58) 

which can be further transformed by using the relation 
between the hypergeometric functions of complementary 
arguments 



F(a, b; c; z) = 
T(c)T(c-a-b) 



F(a,b;a + b- c+ 1; 1 - z) 



T(c-a)T(c-b) 

n ^c-a-6 r ( c ) r (° + & - c ) 

x F(c — a, c — b; c — a — b + 1; 1 — z), 



(59) 



thus resulting in the expression 
I 2 (cos6) + I 2 (-cos6) 



2r(==i) 



-(sin 2 6I) 1 -* 



a + cos ft a — cos ( 

X I T7 + -, 



+ ■ 



\a + cos f 
1 



n 



I a — cos ( 

, 2 



(a + cosft)(a 2 + 2acosft + 1)~ 

— cos 8) 



x F\ 1,— — ;-;sin^7 ] + (cos 



(60) 



Invoking the formula for the derivatives of the hyperge- 
ometric function 



6 



F(a + n, b; c; z) = 



z 1 -* d r ' 
(a) n dz 



oo oo oo 



-[ Z ^F(a,b;c; Z )], (61) x da ds 



dt 



we can rewrite the expression inside the square brackets 
in Eq. (58) as a total derivative. Then the integration 
over a becomes trivial and we finally get the expression 
for the double-source term 

<^ ( % - V) = ^H'^ iL -hx- y\ cos 9] 







x S( P )h [ pt - - e - ia ^ pn + hS(p) 



2 J 



iari\ 



(68) 



r(2fi) 



x \x-y\ z - n {s\n z 6y-- - 



A\x - y 



3—n 



(2-n)(3 - n) 



x F 1 



n — 3 n 



2 ' 2 
which behaves as 



— ; sin 



SG w {x~y) ~ 



4 (^|a;-j/|) 3 
iVvr 2 3 - n 



(62) 



(63) 



After inserting into the integrand the identity J °° dp5(p— 
s — t) = 1 and rescaling the variables s — > sp, i — > one 
can readily perform the integration over s. 

The integration over p results in the table integral 

// \ a / 2 
dxx a - l e- px - q ' x = 2 ( K a {2y/pq), (69) 



thus yielding 
(5G (2) (x-y) 



near n = 3, independent of the cutoff L. 

The second term in Eq. (45) containing one insertion 
of the source current is given by the integral 

oo 

8G {2) {x-y) = j da j d n z S{x - z)nS(z - y) 
o 

x [D(z - y - an) - D(z - x - an)}. (64) 

This expression can be readily computed in the momen- 
tum space where it reads as 



l 



32i r d n p 

^(4tt)™/ 2 J (2tt)"' 



-ip(x-y) 



x(2\ P \y 



J daoM" J dt[t(l - t)]^ [VI - i 



x s\n(atyn)S(v)hv K s- n (a\v\\/ t(l — t)) 



- \p\S(p)Vicos(atpn)Ki^(a\p\y/t(l - t)) . (70) 

The remaining integrals over a are given by the formulas 
6.699.3 and 6.699.4 from the Integral Tables 15 . 
Thus, we arrive at the formula 



VI 



n— 4 



xe lac > n [S(p)nS(p -q)- S(p + q)nS(p)]D{q). (65) 
First, we consider the integral 



x Jdtt—(1 - t) 
o 

5 — n 3 



af \~ - P n -r I 5 ~ 71 
S(p)np.-T^— 



J = j ^ iaqn S{ P -l)D{q) 
N J (2tt)" 



p-q 1 

(P-Cl) 2 ^' 



xF\l, 

xr /3 



->2<-(^H^) 



(66) 



n \ p 1 1 3 - n 1 



t(pn) 2 
Y~ ' 2'~(l-t);p 2 



(71) 



which, upon exponentiating the denominators and inte- 
grating over q, takes the form 



oo oo 



ds 



AV^(47r)"/ 2 J J V~t{s + t) ! i+ 1 





dt 



where the integration over t can be performed by chang- 
ing the variable t = u/(l + u) and comparing the result 
with the integral representation for the hypergeometric 
function 3 F 2 of a certain argument. 

However, one can notice that at n — > 3 the main contri- 
bution stems from the second term in the square brackets 



x e ~P — + ia 7+i P"-4(7+TT 



Thus Eq. (63) can be written as 



5G {2) {x-y) = - 



TVv^Ltt)™/ 2 J (271-)™ 



pt 



d n p 



(67) 



SG {2) (x-y) 
8T (n- 



2r (n-!)r(^) M 3 -" 



x-y 



Nir n +ir (5=2) [(x-y) 2 ]"- 



3—n 



iVvr^rd) 3-n 



-y), 



(72) 



where we restored the dependence on the dimensionful 
parameter u and also used 
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S(x) 



r(n/2) 



2n n/2 ( x 2)«/2- 



(73) 



Combining Eqs. (47), (61), and (70) together, we find the 
overall correction to the fermion propagator 



G(x-y) 



1 + 



32 (fi\x-y\ 



\3-n 



3ir 2 N 



3 - n 



S(x-y), (74) 



from which one can read off the anomalous dimension. 

Remarkably, the latter appears to be still given by 
Eq. (45), as in the case of the original "short-cut" contour 
studied in the previous Section. 

Furthermore, a similar calculation shows that the neg- 
ative anomalous dimension Eq. (45) also pertains to the 
case of the parallel strings which corresponds to choosing 
the upper sign in Eq. (46). 

Taken at their face value, these observations suggest 
that the gauge invariant amplitude Eq. (2) may even be 
largely independent of the choice of the contour T. 



IV. SUMMARY 

In this work, we carried out a direct calculation of the 
previously conjectured form of the physical electron prop- 
agator in such effective QED-like models as the theory of 
the pseudogap phase of the cuprates. In contrast to the 
earlier work, we performed our calculations in the reli- 
able radial gauge and confirmed the result (45) obtained 
in Refs. 10 ' 13 . 

In the course of our analysis, we also investigated the 
dependence of the amplitude (2) on the choice of the 
contour T by considering the case of two (anti)parallcl 
scmi-infinitc strings attached to the end points. The cor- 
responding gauge invariant amplitude is given by Eq. (2) 
with the current Eq. (46) entering the line integral (8). 
Remarkably, the algebraic behavior (1) controlled by the 
same negative anomalous exponent (45) appears to be 
valid for these functions as well. 

In addition to the possible dependence on the choice 
of the contour T (or a lack thereof), the anomalous di- 
mension may strongly depend on the massless fermion 
amplitude in question. For instance, when computed in 
one of the covariant gauges, the amplitude 



G 6 (x-y) 



(0 \ip{x) exp 


iit-Vfdz^A^z) 

y 


4>(y)\o) 


(0|cxp 


X 

V 


|0> 



exhibits a positive anomalous dimension 
16 



3ir 2 N 



(3£ - 2) 



(75) 



(76) 



for any £ > 2/3 13 , including the case of £ = 1 which has 
been claimed 11 to provide an identical representation of 
the original function Go(x) given by Eq. (2). 

However, for any £ ^ the amplitude G^(x) given by 
(75) is not truly gauge- invariant, and, in particular, its 
anomalous dimension computed in a non-covariant gauge 
may differ from Eq. (76). (for an extended discussion of 
this subtle point, see 13 ). For instance, when computed 
in the radial gauge applied in this paper the anomalous 
dimension of the function G^ (x) turns out to be indepen- 
dent of £ and is given by Eq. (45). 

This makes it clearly impossible to substitute any of 
the surrogate amplitudes G^(x) with £ > 2/3 (e.g., 
Gi(x), as in Ref. 11 ) for the original one, Go(x), which 
is the only truly gauge-invariant member of the family of 
functions (75). 

Evidently, the negative anomalous dimension mani- 
fested by the function (2) contradicts the anticipated be- 
havior of a viable candidate to the role of the physical 
electron propagator, since in all of the previously dis- 
cussed effective QED 3 -\ike models the repulsive electron 
interactions are expected to result in further suppression, 
rather than enhancement, of any amplitude describing 
propagation of physical electrons. 

In particular, the algebraic decay of the gauge- 
invariant fermion amplitude (1) would only result in the 
sought-after Luttingcr-likc (stronger-than-linear) vanish- 
ing density of states u{e) ~ |e| 1+?) if r\ were positive. 
By the same token, a pseudogap theory can only recon- 
cile with the experimentally established absence of well- 
defined nodal quasiparticles if a branch-cut singularity of 
the electron propagator that occurs at p 2 ^ = appears to 
be weaker (not stronger!) than a simple pole. 

We defer a further discussion of the construction of the 
physical electron propagator until future work (see, how- 
ever, Refs. 10 ' 13 for an alternate form which demonstrates 
a faster-than-algebraic decay, thus further diminishing 
the chances that the conjecture about the Luttinger-like 
behavior in QED 3 may still be "right, albeit for a wrong 
reason" ) . 

Instead, we suggest that the negative anomalous di- 
mension (45) of the heuristically chosen gauge-invariant 
amplitude (2) may pertain not so much to the physical 
electron propagator per se, but rather to the vertex cor- 
rections which also control the behavior of various gauge- 
invariant two-particle amplitudes ("susceptibilities"). 

In this regard, we quote the earlier result of Ref. 17 
obtained for the susceptibility associated with the four- 
fcrmion scalar vertex 



< 0\i/j(x)ip(x)ip(y)ip(y)\Q >cx 



1 



x-y 



4-(64/37r 2 W) 



(77) 



which features a negative anomalous dimension 2r]. In 
the context of the QED3 theory of the pseudogap phase 
of the cuprates, the formula (77) describes the divergence 
of the staggered spin susceptibility at the antiferromag- 
netic ordering vector Q = (71-, 7r) 18 . 
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We emphasize that one encounters the above problem 
with the unphysical (slower than oc 1/x 2 ) decay of the 
amplitude (2) only in the massless case, while for a finite 
fcrmion mass this function decays as oc e _m l x L 

In the case of the QED 3 theory of the pseudogap 
phase, it has been argued that one may indeed expect 
a dynamical mass generation corresponding to the in- 
trinsic instability towards a spin and/or charge density 
wave ordering 8 ' 11 . 

The question remains, though, as to whether or not the 
chiral symmetry breaking can at all occur for the physical 
number of fermion flavors (TV = 2). Even in the fully 
Lorentz-invariant situation there exist some analytical 19 
and numerical 20 results which suggest the upper bound 
N cr < 2 for the critical number of flavors below which 
the chiral symmetry gets broken. 

In the (non-Lorentz-invariant) QED$ theory of the 
pseudogap phase of Refs. 4 ~ 8 , the role of the strong spatial 
anisotropy of the quasiparticle dispersion and, in particu- 
lar, its effect on a possible universality (or a lack thereof) 
of the critical value of N cr still remain to be ascertained 
(see 21 for a discussion of the weakly anisotropic case). 

It is worth mentioning, however, that in the extreme 
non-Lorentz-invariant limit of the QED^-like theory de- 
scribing the problem of layered graphite the estimated 
value of N cr was found to be even lower than in the orig- 
inal Lorentz-invariant case 22 . 

We conclude by stressing that the problem of con- 
structing the true physical electron propagator in the 
effective massless QED-\ike theories still remains unre- 
solved. Nevertheless, our calculation confirms once and 
for all that the naive ansatz (2) is not up to the job, 
thereby eliminating the current basis for the theoretical 
predictions of the Luttingcr-like behavior in the under- 
doped cuprates 5-7 . 

It is, however, conceivable that, while being unrelated 
to the actual behavior of the electron propagator, the 
negative anomalous dimension (45) manifests the same 
properties of the gauge invariant vertex corrections as 
those exhibited by the physically relevant two-fermion 
amplitudes such as Eq. (77). 
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APPENDIX A: FOCK-SCHWINGER PHOTON 
PROPAGATOR 

In this Appendix we demonstrate that the line inte- 
gral in Eq. (2) vanishes in the so-called radial or Fock- 
Schwinger (FS) gauge. We also derive the expression for 
the photon propagator in this gauge. 

The FS gauge is defined as 

(x-xo^A^x^o) =0. (Al) 

In contrast to such widely used gauges as the Landau 
0M M (a;) = 0, the Coulomb d l A l (x) = 0(i = 1,2) and 
the axial n M A M (x) = ones, the FS gauge may break 
the translational invariance because of the presence of a 
fixed point xq. However, an important advantage of the 
Fock gauge is the explicit relation between the potential 
A f _ l (x, xo) and the field strength F^ v 

l 

A l _ l (x,xo)= J daa(x - xq) v F VIJl {a{x - xq) + xo,xo)- 
o 

(A2) 

In order to derive Eq. (A2) we differentiate Eq. (Al) 

An(x,x ) + (x - x ) l/ d l _ l A v (x,x ) = 0, (A3) 

and then use F^ v = d^Ay — d v A^ 7 to write 

A^(x,x Q ) + (x - x y(F^ u (x,x ) + d v A^(x,x )) = 0. 

(A4) 

Upon changing the variable x — > a(x — x ) + x the last 
expression turns into 

-^-[aAJa(x - x ) + x Q ,x )] 
da 

= a(x - x ) v F v ^{a(x - x ) + x , x ). (A5) 

Integrating over a and using the boundary condition 
4,1(10,10) = (see, Eq. (A3)) which assumes the reg- 
ularity of A^^XjXo) at x = xo, we arrive at Eq. (A2). 
Note that the boundary condition A ll (xo,x ) = is es- 
sential for eliminating a residual gauge freedom which 
remains even after imposing the gauge condition (Al). 

Indeed, in addition to the solution (A2) Eq. (Al) can 
be satisfied by any function 

Al(x,x ) = d£f(x-x ), (A6) 

where / is an arbitrary homogeneous function of 2 — 
xo of zero degree. Any such function would necessarily 
be singular at x = x , though. Hence, the regularity 
condition at x — x can be used to fix the residual gauge 
freedom in (Al). 

Under the translation U~ 1 F lIiV {x)Ua = F fliy (x — a) the 
solution (A2) transforms as 
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u a lA n(x,x )U a = A^(x,x ) = A^(x -a,x - a). (A7) 

When expressed in terms of the center mass X — (x + 
y)/2 and the relative x = x — y coordinates the line inte- 
gral in Eq. (2) takes the following form 



I(x,X;x 



X 

= J dz p 



A»{z) 



(x - yY J daA p (a(x -y)+y, x ) 
o 

1/2 

"■•'•()). 

-1/2 



= x^ J daA li (ax + X,: 



(A8) 



Under translations, Eq. (A8) transforms according to the 
rule: I(x, X; Xq) — I(x, X — a; xq — a). 

We can further restrict the gauge condition (Al) by 
choosing the fixed point x a at the center of mass, i.e., 
xq = X 



{x-XfA ll {x,X) = 



(A9) 



(hereafter, we simplify the notation A p (x,x = X) = 

It can be readily seen that in the gauge (A9) the line 
integral vanishes, i.e., I(x,X;X) = 0. Indeed, from 
Eqs. (A8) and (A2) we obtain 

1/2 

I(x, X; X) = x" J daA^ax + X) 

-1/2 
1/2 1 

^x^x u j da J dPpF^(a/3x + X) = (A10) 

-1/2 

due to the antisymmetry of F vp . 

Furthermore, performing a translation with a = X we 
can cast the gauge condition (A9) in the form 



(x-X)i*A li (x-X) = 0, 
which is identical to 

x»A p {x) = 0. 



(All) 



(A12) 



Next, we derive the photon propagator in the gauge 
(A12) where Eq. (A2) reads as 



A p (x) = J da ax v F vp (ax). 
o 



(A13) 



Thus, we find 



D™(x,y) = (0\TA„(x)A„(y)\0) 



1 

= J dad(3afix° 



a y p {0\TF^(ax)F pv {Py)\0). (A14) 



Since the field strength F pv is a gauge invariant quantity 
the correlator (0\TF ap (x)F pu (y)\0) can be calculated in 
any gauge, including, e.g., the Feynman one where it 
becomes 

x°yf>(0\TF^(x)F pv (y)\0) = x'tf (<W^ - <WW 
-5„ v ffffl + & ap dZd»)D{x-v) 

= H^(x,y)D{x-y). (A15) 

Here D(x) is the photon propagator in the Feynman 
gauge 



D(x) = 



4 m 3 ~'T(^) 



A = 



JV ^ 2 



(A16) 



With the use of the relation x p d p = \x\d\ x \ the operator 
H pu {x,y) can be written in the form 

H pv {x,y) = 5 pv d\ x \d\ y \\x\\y\ - d p x u d\ y \\y\ 

-dy yp d lxl \x\ + d*dyx- y . (ait) 

Now, making use of the identity 
l l 
d\x\ J ' da\x\f(ax) = d\ x \ J da\x\f{a\x\x) 
o o 

\x\ 

= d\ x \J d/3f((3x) = f(x), (A18) 



where x — x/\x\ is the unit vector, we obtain 
D pv (x,y) = {Q\TA p {x)A v (y)\Q) 

dad(3 H pu (ax, /3y)D(ax — (3y) 
l 

= H pl ,{x,y) j dad(3D((ax - (3y) 2 ) = S pv D(x - y) 
o 

l l 

- J da d*x v D(ax - y) - J d/3 dly p D{x - 0y) 

o o 
i 

+ J dad(3 d x p d v u x ■ yD(ax - 0y) (A19) 
o 

(cf. Refs. 23 < 24 ). 

While the first term in (A19) depends solely on x, the 
others have a more complicated dependence, hence the 
photon propagator computed in an arbitrary FS gauge is 
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not necessarily translationally invariant. Moreover, the 
last term in Eq. (A19) displays a divergence at n = 3 
which forces one to use the method of dimensional regu- 
larization when computing various quantities. As pointed 
out in Rcf. 24 (also, see Eq. (40)), the divergence of the 
free FS gauge propagator at n = 3 dimensions is, in fact, 
necessary for obtaining the correct results. 

Finally, in order to return to the gauge (A9) we replace 
x ^ x — X, y ^ y — X m (A19), thus obtaining the FS 
gauge propagator D^(x,y) as a sum of the four terms 
(here the arguments x and y are unrelated to the end 
points in the line integral) 



D^(x,y) = S^D(x-y), 



(A20) 



D${x,y) 



j dad*(x - X) v D(ax - y + (1 - a)X), 
o 

(A21) 



D$(x,y) = -J dpdUy - X)»D{x -(3y-(l- 0)X), 
o 

(A22) 

i 

y) = J dad(3d*dy(x - X) ■ (y - X)D(ax - 0y 



-(a-/3)X). 



(A23) 
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